Fault tolerant quantum computation with very high threshold for loss errors 
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Many proposals for fault tolerant quantum computation (FTQC) suffer detectable loss processes. 
Here we show that topological FTQC schemes, which are known to have high error thresholds, are 
also extremely robust against losses. We demonstrate that these schemes tolerate loss rates up to 
24.9%, determined by bond percolation on a cubic lattice. Our numerical results show that these 
schemes retain good performance when loss and computational errors are simultaneously present. 
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One of the most important achievements of quantum 
information theory was the discovery of FTQC: arbitrar- 
ily precise, scalable quantum computations can be per- 
formed using error-prone components, as long as the error 
rate is below a certain threshold [1, 2 . Current thresh- 
olds, around Pthres ^ 10~^ per elementary operation 
[nilll, come tantalisingly close to the experimental state 
of the art. Of particular interest is the beautiful proposal 
developed by Raussendorf, Harrington and coworkers [4]- 
[7 using ideas from topological quantum computing [8^ 
and with only nearest neighbour gates. 

Many proposals for quantum computing suffer qubit 
loss, such as photon loss, atom or ions escaping from 
traps, or, more generally, the leakage of a qubit out of the 
computational basis in a multi-level system. However, 
such errors can be detected and located without affecting 
the state of the remaining qubits. It is therefore expected 
that appropriately tailored FTQC schemes can tolerate 
higher rates for loss errors than for unlocateable errors 
[9j. Certain schemes have loss thresholds of pioss < 0.5 
p^m [TT] , although it is unclear how they perform in the 
presence of computational errors (unlocated errors acting 
within the computational subspace). FTQC proposals 
tolerant to both error types have thresholds of pioss ^ 
3 X 10-3 and Pcomp < IQ-"^ [lj[13]. 

In this Letter, we describe a FTQC scheme which tol- 
erates both loss and computational errors with very high 
thresholds. The thresholds are characterized by a con- 
tour in (pioss^Pcomp) parameter space passing through the 
points (0.249,0) and (0,0.0063). This represents an im- 
provement of almost 2 orders of magnitude over earlier 
results [12] [14]. Our approach requires only that losses 
are detected at the final readout stage. 

Our scheme combines methods from Raussendorf's 
topological scheme [IHS] and our previous work on loss 
tolerance in surface codes [151 [16]. As in the latter, the 
loss tolerance threshold follows from the bond percolation 
threshold on the relevant lattice, here the cubic lattice 
[17 . To this end, we develop and implement a new clas- 
sical algorithm to analyse the syndrome in the presence 
of both loss and logical errors. Monte carlo simulations 
of the resulting FTQC scheme yield an estimate of the 




FIG. 1. (a) Unit cell showing qubits centered on faces and 
edges of the primal lattice. Heavy lines indicate CPHASE 
gates, (b) A primal identity gate showing logical input oper- 
ators, Xin and Zin, and output operators, Xout and Zout, and 
the correlation surfaces, ax and az, that interpolate between 
them. Also shown are physical qubits (spheres) living on the 
correlation surface between the Z operators, and the defect 
regions, D, embedded in the bulk vacuum, V. (c) Failed 
syndromes for error chains on the primal, E, and dual, E* , 
lattices, and candidate correction chains, C and C*. 



(Pioss^Pcomp) threshold contour. This contour exhibits 
a moderate tradeoff between pioss and Pcomp, indicating 
that the code is robust against both kinds of error. 

We first give an overview of Raussendorf's scheme for 
FTQC. The scheme is based on the one-way quantum 
computer: that is, a cluster state |C)/:, of many physical 
qubits located on the faces and edges of a cubic primal 
lattice £, is prepared by first preparing every physical 
qubit in the state |+), and then applying CPHASE gates 
between neighbouring qubits on the lattice [see Fig.[TJa)]. 
It is convenient to introduce the dual lattice, £*, whose 
vertices, edges, faces, and cubes correspond, respectively, 
to the cubes, faces, edges, and vertices of C. Physical 
qubits reside only on the faces and edges of these lattices. 
Computation proceeds by making a sequence of single 
qubit measurements on |C)/:, aided by classical process- 
ing of the measurement outcomes. Correlations between 
the measurement outcomes give rise to fault tolerance. 

C is subdivided into three regions, denoted F, D and 
S. Qubits in V and D are measured in the X and Z 
bases, respectively. Qubits in S are measured in either 



the Y OT {X -\- Y)/^/2 basis. D comprises a collection of 
thick structures with some characteristic transverse di- 
ameter {'defects'), embedded in V [Fig. [ijb)]. Each pair 
of defects encodes one logical qubit. As the D qubits are 
measured in the Z basis, one can omit them altogether, 
which we assume for the remainder of this paper. 

The topology of the braiding of defect regions effect 
certain Clifford gates between the logical qubits, whilst 
measurement results in V provide topologically-protected 
fault tolerant error correction of the logical Clifford gates 
to arbitrarily high accuracy. S comprises a collection of 
well-separated single qubits, spread out among the de- 
fects. The S qubits are used to introduce noisy encoded 
magic states into the circuit, which can then be distilled 
using the very accurate logical Clifford gates in V and 
D, implementing a universal gate set [18 . 

Magic state distillation tolerates a rather large amount 
of noise, and so it is found that the thresholds for both 
loss and computational errors are set by the correspond- 
ing thresholds in the bulk topological region, V and D. 
We therefore discuss error correction in V and D in more 
detail, focusing on the implementation of a fault tolerant 
identity gate. The reader is referred to [HllQ] for details. 

The layout of the identity gate is show in Fig. [ijb). 
The boundaries of D are positioned such that faces of 
C lie just inside the defect regions. The alignment of 
the defects specify a 'simulated time' axis, such that the 
identity gate maps the logical qubit from the input (rear) 
plane, /, to the output (front) plane, separated by an 
integer number of unit cells. Within these planes, the 
qubit is encoded in a surface code [8 , containing two 
holes, which coincide with the intersection of the defect 
regions with the corresponding plane. A single logical 
qubit is associated to the pair of holes, with encoded 
logical operators X and Z as indicated on Fig. [ijb). 

To demonstrate the operation of the noiseless gate, 
it suffices to show that the input logical operators Xin, 
Zin are mapped to the output operators Xout, ^out, un- 
der the action of the CPHASE gates and single qubit 
measurements [20]. To this end, we introduce the clus- 
ter stabilizer operators, Kf = Xf ^^^Qr Z^, associated 
with each face / of £ or £*, where df denotes the 
qubits at the edges of /. The initial state of the clus- 
ter satisfies Kf\C)jc = \C)jC. Consider the surface a^ 
illustrated in Fig. [ijb), consisting of faces in C The 
product of the cluster stabilisers on a 2 yields i^(cr^) = 
(^ p^^_ Xf (2)eGacr- ^e7 where da 2 denotes the qubits at 
the boundary of a^. Since (^ 
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follows that measurements in the bulk effects the opera- 
tor mapping Zin -^ ±^in ^ ^out, where the sign is given 
by the parity of the measurement outcomes Uj^fj-Xf. 

Similarly, the product of cluster stabilisers act- 
ing on the surface aj^, consisting of faces in £*, 
yields K{ax) = ®f*ea^ ^r ®e*eda^ ^e*^ and since 




(8)e* 



e*Gc^cr. 



Xp* — Xu 



) Xoiit it follows that measurements 



FIG. 2. (a) A correlation surface suffering loss is repaired by 
deforming around lost qubits (white spheres), by multiplying 
by the corresponding check operator. An affected logical in- 
put or output operator is similarly deformed, (b) In the bulk, 
super-check operators are formed by products of elementary 
check operators that are affected by a loss. Error chains ter- 
minate within a super-check operator. 



on ax effect the mapping Xin -^ i-?in '^ A'out • Finally, 
measuring the qubits in / leads to the desired trans- 
formation (up to an unimportant Pauli frame update) 
Zin -^ i^out and Xin -^ ±Xout- 

Error syndromes are revealed by correlations of mea- 
surement outcomes. In particular, considering the prod- 
uct of i^/'s centred on the faces of a unit cube c of C 
leads to the parity check Pc = Uf^QcXf. Absent any er- 
rors, the parity checks satisfy Pc = +1. Failed parity 
checks, i.e. cubes which instead satisfy Pc = —1, reveal 
the end points of chains, E, of Z errors (see Fig. [lie)). 
These chains reside on edges of £* , and are corrected by 
inverting the recorded measurement outcomes on some 
correction chain C (i.e. Xf -^ —Xf for f e C) sharing 
the same endpoints as E (i.e. dC = dE). This correction 
process works provided the combined chain E-\-C neither 
winds around a defect, or joins two defects together. 

Finding a suitable correction chain C requires an algo- 
rithm for pairing failed parity checks, then finding a suit- 
able path in the lattice between each pair of syndromes. 
Choosing C optimally is computationally difficult, since 
it is equivalent to minimising the free-energy of a partic- 
ular spin-glass model [6l \2V . However efficient heuristic 
methods for decoding error syndromes have been devel- 
oped, including those using Edmonds' perfect matching 
algorithm [2J^ 22, , and belief propagation [24 . 

Realistic error models depend on the particular im- 
plementation under consideration. As in [7] we assume 
single qubit depolarising errors take place during prepa- 
ration, storage, and measurement, with rates pp,ps and 
Pm respectively, whilst CPHASE gates suffer two-qubit 
depolarising errors with rate p2- For numerical calcula- 
tions, we assume these errors are iid and occur with equal 
probability, Pcomp, i.e. Pp = Ps = Pm = P2 = Pcomp- 
These parameters may vary, with modest changes in the 



threshold [7]. For this error model, previous results es- 
tablished an error threshold of Pcomp < Pt = 0.0058 
[7]. Note that without losses, exploiting error correla- 
tions within each sub-lattice improves the threshold to 
Pt = 0.0075 ^. A further improvement may be possible 
by exploiting correlations between C and £* [24 . For 
simplicity, we do not exploit correlations, at the cost of 
a marginal reduction in the threshold. 

We assume losses are iid with total rate pioss, and only 
occur either before or after all CPHASE gates have acted 
upon each qubit, but not at intermediate times. This is 
consistent with e.g. optical lattices where loading losses 
[25] are much higher than storage losses [26 , or with a 
photonic implementation where source and detector in- 
neficiencies are expected to dominate. CPHASE gates 
acting on a lost input are assumed to implement the iden- 
tity operation, plus depolarizing noise at rate p2^ on the 
remaining qubit. Losses need only be detected at the 
measurement step, so loss detection is permitted to be 
a destructive process. Loss detection may be imperfect, 
but such errors are accounted for in pm or pioss • 

To deal with these losses, we adapt the Edmonds' 
matching approach, paralleling [151 116]. Suppose qubit 
q is lost. The two parity checks Pc^ and Pc^ on the ad- 
jacent cubes Cq and Cqf are incomplete so cannot detect 
error syndromes. Likewise any correlation surface a that 
depends on q is damaged so cannot mediate logical gates. 

We recover from this damage in two ways. Firstly, 
multiplying K{(j) by K{dcq) = ^f^Q^ ^fi yields a new 
surface d that avoids q [see Fig. ^a)], i.e. K{d) = 
K{a)K{dCq) is independent of q. Then a determines an 
equivalent mapping of input operators to output opera- 
tors. This procedure can be iterated for all of the lost 
qubit s, and succeeds as long as there is some 2D corre- 
lation surface that percolates between the logical qubit 
operators whilst completely avoiding lost qubit s. In the 
limit where the defect dimensions and spacings between 
defects becomes large, such a surface can be found pro- 
viding the loss rate is lower than the cubic-lattice bond 
percolation threshold, which is approximately 0.249 [17]. 

Secondly, the product of the two damaged check oper- 
ators Pq = Pc Pc' is independent of g, so the sign of Pq 
yields information about dE^ as illustrated in Fig. ^b). 
Thus for each instance of losses, we form new 5'^per-check 
operators that are insensitive to the losses. The lattice 
of check operators is then no longer cubic, so we perform 
the syndrome matching on an irregular lattice. 

In the context of surface codes, Edmonds' match- 
ing has typically been used to find maximum likeli- 
hood correction chains, C, that maximise the probability 
P{C\dE) [4l[22]. It is simple to show that such chains are 
also of minimum length [22] . For a typical pattern of loss, 
it turns out that there are very many minimum length 
chains. Furthermore, as discussed in [15] some minimum- 
distance pairings have many more matching chains than 
other minimum-distance pairings, and so they are a pri- 



ori much more likely. Accounting for this degeneracy in 
the pairing is computationally cheap, and leads to modest 
improvements in the threshold. In the 3D FTQC scheme 
on an irregular lattice, the degeneracy of matchings may 
also be computed efficiently by modifying Dijkstra's al- 
gorithm for finding shortest paths on a graph. 

To establish quantitative error thresholds in the pres- 
ence of losses, we perform Monte Carlo simulations on 
L^ cubic lattices with periodic boundary conditions. For 
each simulation we do the following: 

1. Generate an instance of errors E and losses. 

2. From losses, derive P and d. If losses percolate the 
cubic lattice then failure (no such d exists). 

3. From E^ derive dE on the irregular lattice. 

4. From dE^ compute the correction chain C. 

5. Compute the homology class 1-Le+c (count the 
number of intersections oi E ^ C with <j, mod 2). 

6. If 1-Le+c = then SUCCESS. Otherwise failure. 

At a given loss rate, we generate 198, 000 instances 
of errors whilst varying Pcomp and L = 6,8,... 16. For 
each loss rate, we fit a Taylor's expansion of a univer- 
sal scaling function Pfaii(pioss,Pcomp, L) ^ a ^ bx ^ cx'^, 
where x = (pcomp —Pt)L^^^j with fitting parameters pt 
(the loss-dependent computational error threshold), u 
(the scaling exponent), a, b and c. 

Fig.|3]is the central result in this paper, and shows pt at 
different values of pioss (i"ed circles). The shaded region 
indicates the correctable region of parameter space; if 
Pcomp < Pt^ then errors are almost always correctable (in 
the limit L -^ oo). The solid curve is simply a quadratic 
fit to the computed values of pt for < pioss < 0.15 [27] . 
and extrapolates to pioss = 0.252 ± 0.005 at Pcomp = 0, 
which agrees with the percolation threshold Pperc = 0.249 
[17] (green diamond) on a 3D cubic lattice. 

The loss-tolerant thresholds above rely on properties 
of the FTQC scheme in the bulk V and D regions, well 
away from S qubits. Special attention must be given to 
the effect of losses affecting S qubits (and also to nearby 
V and D regions), which are used to inject encoded magic 
states into the circuit. Computational errors close to 
the S qubits increases the effective noise on the encoded 
magic states to ps ~ Qpcomp^ but the threshold for magic 
state purification is so high that the overall threshold 
is still set by that of the bulk [4]. This remains true 
with losses, and we now describe a scalable post- selection 
method which, whilst profligate, demonstrates that losses 
near S qubits do not limit the thresholds. 

If a loss occurs either on or near an S qubit, we dis- 
card the corresponding encoded magic state by mea- 
suring the encoded state, i.e. fusing the neighbouring 
defect strands together. This may be performed in 
the bulk region, provided the error rates are below the 
corresponding bulk thresholds. Thus, we require that 
the only magic state qubits which are injected into the 
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FIG. 3. Phase diagram showing correctable region of param- 
eter space. Red circles are numerically calculated computa- 
tional error thresholds, pt at different loss rates, for an error 
model in which pp — ps — pm — P2 — Pcomp. The green 
diamond is the 3D bond percolation threshold 17 . The blue 
square is the threshold computed in 7 , which ignores degen- 
eracy in the matching algorithm. 



circuit are those for which no losses occur within a 
d X d X d volume centred on the S qubit itself. Then 
the post-selected error rate on the encoded magic states 
is PS ^ 6pcomp +Pfaii(pioss,Pcomp,c^). Bclow the bulk er- 
ror threshold, pfaii is exponentially suppressed in d. The 
magic state distillation threshold p^ < 1/a/35 [18 fixes d 
to be a (small) constant and so the additional overhead 
for this post-selected scheme, ^ (1 — pioss)^ , is indepen- 
dent of the size of the algorithm, ensuring scalability. 

In this letter we have established high thresholds 
FTQC for both computational errors and loss er- 
rors. This result, which sets an important bench- 
mark for experimental implementations, follows from the 
large redundancy of the surface codes that underpin 
Raussendorf's FTQC scheme. This redundancy serves 
two complementary purposes: to enable deformed corre- 
lation surfaces to mediate logical gates between encoded 
qubits, and to enable expanded parity check operators to 
identify end points of computational error chains. 

A number of lines of enquiry follow from this work. 
We assumed deterministic CPHASE gates, so relaxing 
this to allow heralded, non-deterministic two-qubit gates 
in the cluster preparation [28^30] is open (but see [31 ). 
Also, we expect that losses which occur during cluster 
construction will reduce thresholds further, since these 
induce additional Z errors on neighbouring qubits. Since 
these are located errors, percolation phenomena still de- 
termine the loss threshold, so we anticipate this will re- 
main substantially higher than the computational error 
threshold. Finally, it is important to develop methods 
to deal with losses on S qubits with lower overhead than 



that presented here, perhaps by dynamically routing de- 
fect regions and S qubits around any revealed losses. 
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